Interface resistance of disordered magnetic multilayers 
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We study the effect of interface disorder on the spin- 
dependent interface resistances of Co/Cu, Fe/Cr and Au/Ag 
multilayers using a newly developed method for calculating- 
transmission matrices from first-principles. The efficient im- 
plementation using tight-binding linear-muffin-tin orbitals al- 
lows us to model interface disorder using large lateral super- 
cells whereby specular and diffuse scattering are treated on 
an equal footing. Without introducing any free parameters, 
quantitative agreement with experiment is obtained. We pre- 
dict that disorder reduces the majority-spin interface resis- 
tance of Fe/Cr(100) multilayers by a factor 3. 



When two layers of magnetic material are separated by 
a non-magnetic spacer layer, the electrical resistance of 
the system depends strongly on whether the magnetiza- 
tion directions are aligned parallel or anti-parallel. This 
effect is known as giant magnetoresistance (GMR) jjj. 
The huge interest [§-f|] in the physics of GMR is largely 
driven by the wide application potential of the effect, 
which has already been realized in magnetic recording 
heads. 

GMR can be observed in a number of different mea- 
suring configurations. The current-in-plane (CIP) con- 
figuration is experimentally the simplest and is what is 
used at present in applications. However, for gaining 
a better understanding of the underlying physics, the 
current-perpendicular-to-the-plane (CPP) configuration 
H|5|-|| is preferred because of its higher symmetry which 
should make it easier to understand and because of higher 
MR ratios. 

The factors usually considered in theoretical treat- 
ments of GMR are the potential steps encountered by 
electrons passing from one material to another, impurity 
scattering in the bulk of the layers, and defect scatter- 
ing at the interfaces |2|-||]. There has been a great deal 
of discussion about the relative importance of these in- 
gredients and their spin dependence, which cannot be 
resolved solely on the basis of model calculations which 
include these effects in parameterized form. Once the 
question has been suitably posed, however, detailed elec- 
tronic structure calculations can be used to resolve the 
issue quantitatively. For example, the effect of potential 



steps and their microscopic origin could be established in 
this way 

In this Letter we wish to address the relative role of 
specular and diffuse interface scattering. This has been 
studied by a large number of authors but so far only using 
simple models which do not allow for detailed quantita- 
tive analysis of specific materials We focus on 
the interface resistance of the resistor model which de- 
scribes the observed thickness and layer dependence of 
CPP-GMR remarkably well ||. Because it turns out to 
be strongly spin-dependent and dominates the magne- 
toresistance for layer thicknesses which are not too large, 
the key to understanding CPP magnetoresistance lies in 
understanding the origin of the interface resistance. The 
methodology which we have developed allows us to in- 
clude specular and diffuse scattering on an equal footing 
without introducing any arbitrary fitting parameters. 

Explicit expressions for the interface resistance were 
derived by Schep et al. |l5| ] in terms of the transmission 
matrix T which describes how the electronic structure 
mismatch at an A/B interface affects electron transport. 
In the limit in which there is no coherent scattering be- 
tween adjacent interfaces, presumably due to sufficiently 
strong bulk scattering, the interface resistance is given 

by 
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where T^ v are the probabilities for eigenstate /i in mate- 
rial A to be transmitted through the interface into the 
eigenstate v in material B, the sum is over the Fermi 
surface, and e 2 /h Na(b) is the Sharvin conductance of 
material A(B). In Ref. |[3| transmission matrices and 
interface resistances were obtained for ideal Co/Cu in- 
terfaces using a first-principles FLAPW based embed- 
ding technique. These, and similar results obtained by 
Stiles and Penn H], demonstrated that a combination 
of spin-independent bulk scattering and strongly spin- 
dependent specular interface scattering arising from the 
spin-dependence of the band mismatch can account for 
the observed spin-dependence of the interface resistances. 
These results are at odds with the common wisdom that 
metallic heterointerfaces cannot be perfect due to un- 
avoidable roughness and/or interface alloying. Indeed, 
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for the one case of Co/Cu(lll) interfaces for which direct 
comparison could be made with experiment, the agree- 
ment though reasonable, was not perfect. We therefore 
address the following questions: Why does a calculation 
for a perfect interface agrees so well as it does for a sam- 
ple produced by sputtering? Can theory and experiment 
be brought to even better agreement by taking into ac- 
count disorder? Is the finding that specular interfaces 
are a reasonable first order approximation generic or a 
coincidence found only for the Co/Cu system? 

The FLAPW-based method used in Ref. |l3| was com- 
putationally too demanding to allow interface disorder to 
be treated. Starting instead with the more efficient sur- 
face Green's function method [flEf implemented with a 
tight-binding linear muffin tin orbital basis jL6), we can 
now calculate the transmission and reflection matrices 
needed in the Landaucr-Buttikcr formulation of trans- 
port theory ]l7j , but now for much larger systems. In this 
Letter, we present the results of calculations for Co/Cu, 
Fe/Cr and Ag/Au layered systems in which we model in- 
terface disorder by means of large lateral supercells. The 
electronic structure is determined self-consistently within 
the local spin density approximation. To model the in- 
terface, we randomly distribute the appropriate concen- 
tration of different atoms within lateral supercells |l^] 
containing as many as 10 x 10 atoms. For the disordered 
layers the potentials are determined self-consistently us- 
ing the layer CPA approximation |l5| ]. The calculations 
are carried out with a k|| mesh density equivalent to 3600 
k|| mesh points in the two dimensional Brillouin zone 
(BZ) of a 1 x 1 interface unit cell. The numerical error 
bar resulting from this sampling is smaller than 0.2% of 
the conductance. The interface resistances calculated for 
Co/Cu (100) and (111) in the clean limit using Eq. 
agree with those obtained by Schep using an entirely dif- 
ferent code to within about 0.1 ffhn 2 or 5%. 

In the presence of defects, the conductance can be ex- 
pressed as the sum of a ballistic part and a diffuse part; 
the transmission matrix elements between two Bloch 
states with the same k|| correspond to ballistic scatter- 
ing, those between two Bloch states with different k|| to 
diffuse scattering. 



TABLE I. Results of calculations 



system 


roughness 


R maj (ffim 2 ) 


R min (ffim 2 ) 


Au/Ag(lll) 


clean 


0.094 


0.094 


Au/Ag(lll) 


2 layers 50-50 alloy 


0.118 


0.118 


Au/Ag(lll) 


exp. 


0.100 ±0.008 


0.100 ±0.008 


Co/Cu(100) 


clean 


0.33 


1.79 


Co fecp /Cu(lll) 


clean 


0.60 


2.24 


Co/Cu(lll) 


clean 


0.39 


1.46 


Co/Cu(lll) 


2 layers 50-50 alloy 


0.41 


1.82±0.03 


Co/Cu(lll) 


exp. 


0.26±0.06 


1.84±0.14 


Fe/Cr(100) 


clean 


2.82 


0.50 


Fe/Cr(100) 


2 layers 50-50 alloy 


0.99 


0.50 



The calculated results are shown in Table |. The 
Au/Ag interface has fcc(lll) texture and the interface 
roughness is estimated to be at least two layers thick 
in the MSU samples |Q. This makes it very suit- 
able for testing our method without the complicating 
factor of spin dependence. The interface resistance we 
find for a clean Au/Ag interface based on Eq. ([!]) is 
0.047 ffim 2 which is very close to the experimental value, 
0.050 ± 0.004 ffim 2 . The resistance of a Au/Ag interface 
becomes 0.059 ffhn 2 when the interface contains two lay- 
ers of Auo.sAgo.s alloy. The uncertainty arising from us- 
ing different alloy configurations within the 6x6 unit cell 
is less than 1%. 

We calculate the interface resistance of Co/Cu inter- 
faces for both (100) and (111) orientations. The lattice 
constant used for fee Co/Cu is 3.549 A. We focus mainly 
on the (111) orientation as this is the structure which is 
predominantly seen in the experimental samples. The in- 
terface alloy is again at least two atomic layers thick ||. 
We treat the interface disorder as two layers of CoCu al- 
loy modelled using an 8 x 8 lateral supercell. The largest 
uncertainty between different configurations of two layers 
of 50-50 alloy is about 2.5%, which is much smaller than 
the experimental error bar. For interface alloy composi- 
tions ranging from 50-50 to 44-56 the interface resistance 
does not change within our numerical accuracy. With two 
layers of interface alloy, the calculated transmission prob- 
ability for the minority spin electrons decreases by about 
10% bringing the calculated interface resistance into near 
perfect agreement with experiment. We find that disor- 
der gives rise to mainly forward scattering of the electrons 
so that the decrease of the ballistic component is almost 
cancelled by the increase of the diffusive part. This is the 
reason why the calculations for the defect-free Co/Cu in- 
terface 1^,0 were in reasonable agreement with exper- 
iment. The strong diffuse scattering also explains why 
the two-channel resistor model performs so well down to 
relatively thin layers in which bulk scattering should not 
be important. 

The resistance of a Cu/Co interface calculated with 
two layers of 50% interface alloy is 0.41 film 2 for the ma- 
jority spin and 1.82 ± 0.03 ff2m 2 for the minority spin. 
The error bar for the minority-spin results from using a 
finite lateral supercell for modelling disorder and configu- 
ration averaging. The majority spin bands of Cu and Co, 
like the bands of Au and Ag, are so well matched that 
interface disorder has very little effect on the interface 
resistance. We observe (Table I) that there is near per- 
fect agreement for the minority spin (certainly within the 
overall uncertainty of the calculation) but that the calcu- 
lated resistance for the majority spin case which was al- 
ready too large in the absence of disorder is even slightly 
increased by disorder. 



2 




2 4 
Number of Interfaces 



sure it. 



FIG. 1. Differential interface resistance as the number of 
interfaces increase for a disordered Cu/Co multilayer embed- 
ded between Cu leads. 
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With two layers of interface alloy, almost 80% of 
the minority spin current results from diffuse scatter- 
ing. We can see this in a different way by calculating 
the conductance for Cu/Co/Cu, Cu/Co/Cu/Co/Cu, and 
Cu/Co/Cu/Co/Cu/Co/Cu in a manner quite analogous 
to that described above for the single Cu/Co interface. 
To be able to perform this large calculation we had to use 
a smaller supercell (6 x 6) so that the error bar is larger 
than for the 8x8 calculation. In these calculations the 
boundary Cu layers are semiinfinite "leads" and the in- 
terface disorder is two layers of 50-50 alloy. Results are 
insensitive to the individial layer thicknesses, chosen here 
to be 10 atomic layers. In Fig. [jjthe differential resistance 
per interface is shown for all four systems. For the mi- 
nority spins, the deviation from the result obtained for a 
single interface is quite small. The strong diffuse interface 
scattering destroys the phase coherent scattering between 
subsequent interfaces so that Ohm's law holds when the 
number of interfaces is increased. For the majority spin, 
however, the interface resistance does not obey Ohm's 
law but decreases as the number of interfaces increases. 
It appears to saturate at a value of 0.07 ffim 2 which is 
only a third of the experimental value, but consistent 
with Mathon's calculation for a multilayer with random 
layer thicknesses (2(]]. The majority spin potentials are 
so similar that the scattering from a double alloy layer 
is insufficient to break the coherence which is consider- 
ably longer range than for the strongly scattered minority 
spins. We would have to assume that the majority spin 
electrons remain coherent for transport through 4 inter- 
faces in order to obtain a value of the interface resistance 
close to the experimental value of 0.26 ffim 2 . Compared 
to real samples with bulk defects and lateral variations 
in the layer thicknesses, it appears that we overestimate 
the coherence length in the majority spin case [Q. 

For the (111) orientation we also considered an in- 
terface between hep Co and fee Cu. For a clean 



Fe(1.95) 



Co 



(0001) 
hep 



/Cu^ 11 ^ interface both majority and minority 



spin resistances are substantially larger than for the fee 
case and even larger than the experimental values (Ta- 
ble |). In view of the substantial difference predicted 
between the two, it should be interesting to try and mea- 
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FIG. 2. Number of propagating channels in first Brillouin 
Zone, (a), (b) and (c) are for the majority spin of a clean 
Fe/Cr interface, bulk Fe, and bulk Cr,(d), (e) and (f) are for 
minority spin electrons of a clean Co/Cu interface, bulk Co, 
and bulk Cu, respectively. The total number of propagating 
channels per unit cell is given in brackets on top of the fig- 
ures. The grayscale interpolates the number of propagating 
channels per feii -point between zero (white) and four (black). 

The Fe/Cr (non- magnetic) interface calculation is per- 
formed for the bcc(100) orientation, which is the low in- 
dex orientation with the largest spin-asymmetry |]l4| . We 
used a lattice constant of 2.87 A. To model the interface 
roughness a 6 x 6 lateral supercell was used. The uncer- 
tainty from configuration averaging is less than 10%. 

Whereas for Co/Cu the majority-spin band structures 
were well matched, for Fe/Cr the situation is reversed 
and it is the minority-spin electronic structures which 
match well. Using Eq. ([!]), the interface resistance is 
2.82 ffim 2 and 0.50 ffim 2 for majority and minority spin, 
respectively, so that the Fe/Cr interface has a negative 
spin-asymmetry, opposite to that of bulk Fe. The ef- 
fect of disorder is to suppress the interface asymmetry 
rather than enhance it. As was the case for the Co/Cu 
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majority spins, interface disorder has only a small ef- 
fect on the well-matched minority spin channel. For the 
majority spin channel however, the transmission proba- 
bility for a clean interface is very low due to a large band 
mismatch. For a disordered interface, the ballistic contri- 
bution to the conductance can only decrease by a small 
amount but the diffuse component increases enormously 
leading to a large net increase in the transmission. 3% 
Fe impurities in the first Cr layer (or 3% Cr in the first 
Fe layer) increase the transmission probability by more 
than 10%. Two interdiffused atom layers suppress the 
spin asymmetry and the MR efficiently - the interface re- 
sistances resulting from two 50-50 interface alloy layers 
are 0.99 ff2m 2 and 0.50 ff2m 2 for majority and minor- 
ity spin respectively. Vice versa, the interface quality 
is much more critical for a large CPP-MR in the Fe/Cr 
than in the Co/Cu system. 

The qualitative difference between Fe/Cr and Co/Cu 
can be understood using Figs. ||(a-c) and Figs. g(d- 
f) where we show as a function of k|| the number of 
majority-spin propagating channels for the Fe/Cr inter- 
face, bulk Fe, and bulk Cr, respectively, and the num- 
ber of minority-spin propagating channels for the Co/Cu 
interface, bulk Co, and bulk Cu, all in the first BZ. For 
the Co/Cu (111) interface minority-spin states, the Fermi 
surfaces of both Co and Cu occupy a large part of the 
2D BZ so that there are a large number of states with 
the same k|| in both materials which can, in principle, 
propagate in the absence of disorder. Summing over all 
ky , however, the transmission probability of states (com- 
ing from Cu) is only about 60%; the character of the 
bulk states on either side is such that they match poorly. 
Defect scattering tends to reduce the transmission prob- 
ability and thus increases the interface resistance of the 
Co/Cu minority spin channel. On the other hand, we 
can identify two mechanisms by which interface disorder 
decreases the interface resistance in the Fe/Cr majority- 
spin channel by a factor 3. Majority spin electrons with 
small ky are almost completely reflected at the Fe/Cr 
specular interface because the electronic states on both 
sides of the interface do not match well. Defect scatter- 
ing is found to increase the transmission of these electrons 
strongly. Furthermore, for k|| outside of this central area, 
there are no propagating states on the Cr side. Propagat- 
ing modes in Fe with larger k|| , which are totally reflected 
at the specular interface, can be scattered diffusely into 
the center of the BZ where many states are available in 
Cr. 

In summary, we have studied the interface resistance of 
Co/Cu, Fe/Cr and Au/Ag interfaces. Depending on the 
system, interface disorder can increase or decrease the 
interface resistance. For some interfaces such as Au/Ag 
and Co/Cu, the band mismatch at an interface is respon- 
sible for most of the interface resistance. For other sys- 
tems such as Fe/Cr, the interface resistance can dramat- 
ically depend on the interface perfection. For Fe/Cr in- 



terface, the majority-spin interface resistance is reduced 
as much as 70% by interface disorder. Interface disorder 
enhances the spin asymmetry in the Co/Cu system but 
decreases it for Fe/Cr. In the systems considered, the 
diffuse scattering arising from interface disorder breaks 
the phase coherence in high resistance spin channels, but 
not necessarily for the low resistance spin channels. 
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